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ABSTRACT 

We derive a non-parametric CO luminosity function using a FIR and an optical i?-band selected 
sample of the galaxies included in the FCRAO Extragalactic CO Survey. The FIR selected sample 
is defined using the IRAS Bright Galaxy Surveys (BGS; IRAS 60 micron flux density > 5.24 Jy). 
Although our CO sample is not complete, the normalization using the BGS reproduces the IRAS 60 
micron luminosity function in excellent agreement with those found in the literature. Similarly, a B- 
band selected sample defined using the Revised Shapley-Ames (RSA) catalog is used to derive a CO 
luminosity function for a comparison. A Schechter function describes the both derived CO luminosity 
functions reasonably well. Adopting the standard CO-to-H2 conversion factor, we derive a molecular gas 
density of = (3.1 ± 1.2) x IO^HMq Mpc~^ for the local volume. Combining with the measurements 
of the local HI mass density and the helium contribution, we estimate that the total mass density of cold 
neutral gas in the local universe is ilgas — (4.3 ±1.1) x which is about 20% of the total stellar 

mass density fi*. 

Subject headings: galaxies: luminosity function - galaxies: ISM - galaxies: surveys - ISM: molecules - 
radio lines: ISM - cosmology: cosmological parameters 



1. INTRODUCTION 

The star formation history of the Universe is closely 
linked with the evolution of the gas content of the Uni- 
verse. Observations of the distribution and total gas con- 
tent in galaxies and in intergalactic clouds offer some of 
the most important observational constraints for the cos- 
mology and galaxy evolution models. Because stars form 
out of cold, dense gas clouds, the history of galaxy for- 
mation and evolution is also the history of gas accretion 
and conversion into stars. Here we derive the total cold 
gas density for the local volume by deriving the local CO 
luminosity function and combining with the existing esti- 
mates of the neutral atomic gas density. 

A variety of luminosity functions (LFs) and mass func- 
tions are found in the literature: optical (e.g. Marzke, 
Huchra, & Geller 1994), HI (e.g. Zwaan et al. 1997; 
Schneider et. al 1998), and infrared (e.g. Soifer et al. 
1987; Yun et al. 2001). While several groups have investi- 
gated the local neutral hydrogen density, little information 
is available on the molecular gas content in galaxies be- 
cause of the lack of appropriate data. In her pioneer work, 
Verter (1987) derived averaged values of CO luminosity 
and the CO /HI flux ratio based on a maximum likelihood 
probability distribution for 40 galaxies and 47 upper lim- 
its obtained by combining all available data at the time. 
Since the large fraction of the CO measurements used were 
non-detections (upper limits), the derivation of the CO lu- 
minosity function was highly problematic. Nevertheless, a 
general trend of increasing number of CO emitting galax- 
ies per unit volume with decreasing CO luminosity was 
suggested by this analysis. 

Nearly 15 years later, the availability of the extragalactic 
CO data has improved greatly. The FCRAO (Five College 
Radio Astronomy Observatory) Extragalactic CO Survey 



(Young et al. 1995) represents a particularly rich database 
for the investigation of CO luminosity function and molec- 
ular gas content in galaxies. We have constructed a large, 
statistically significant sample of far-infrared (FIR) and 
optical i?-band selected galaxies from this survey and de- 
rived a non-parametric CO luminosity function. Our sam- 
ple galaxies range over 4 orders of magnitudes in CO lumi- 
nosity, and only a small fraction are non-detections. Since 
CO is a tracer of hydrogen molecules, molecular hydrogen 
mass can be derived from CO luminosity. By integrating 
the resulting molecular gas mass function, we then derive 
the molecular gas mass density and the total cold gas mass 
density of the local universe. 



2. SAMPLE SELECTION AND PROPERTIES 

The FCRAO Extragalactic CO Survey ("the Survey", 
herein) is the largest CO survey of galaxies available in 
the literature, containing observations of 300 galaxies at 
1421 positions (Young et al. 1995). Since the Survey 
covers a broad range of Hubble types and contains a large 
number of detections (236 or 79%), it is well suited for a 
statistical analysis and for constructing a CO luminosity 
function. The majority of the survey galaxies were selected 
from the Second Reference Catalog (RC2; de Vaucouleurs 
et. al 1976) or the IRAS database. Most of the galaxies 
in the Survey are spiral or irregular galaxies at declina- 
tion north of —25° and satisfy at least one of the following 
three criteria: < 13, Seo^rn > 5 Jy, or Sioofim > 10 Jy 
(see Young et al. 1995, for a detailed description). The 
majority of the detected galaxies in the Survey were ob- 
served at multiple positions along the major axis, and the 
total CO fluxes were derived by modeling the underlying 
CO distribution. 
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The Survey is not complete in terms of any of the sam- 
ple selection criteria, and the full sample cannot be used 
directly for a statistical analysis in a straightforward way. 
The analysis of the velocity integrated CO luminosity is 
further complicated by the fact that the CO line width 
varies from one galaxy to another (and even one posi- 
tion to another in a given galaxy) while the sensitivity 
achieved also varies depending on the observing condition 
even if identical integration times are used (typically 2-4 
hours). Therefore, it is difficult to address tlic^ sensitiv- 
ity and completeness limit of the Survey in terms of the 
observed CO flux. Instead, taking advantage of the well- 
known tight correlation between the FIR and CO luminos- 
ity (see Young et al. 1995), we utilize the IRAS 60 /xm 
flux density of the individual sample galaxies to define a 
complete sample for a statistical analysis and the deriva- 
tion of the CO LF. We further test the robustness of the 
technique and the selection bias by deriving the CO lu- 
minosity function using the optical i?-band selection (sec 
§ 3.3). Previously Briggs & Rao (1993) successfully de- 
rived the HI mass function for the field galaxies using the 
optical selection, and the radio luminosity function has 
been successfully derived using the FIR selection function 
by Yun et al. (2001). 

Taking advantage of the well defined and well studied 
IRAS Bright Galaxy Samples (BGSs; Soifer et al. 1989; 
Sanders et al. 1995), we adopt the sample selection cri- 
terion of 60 /im flux density limit greater than 5.24 .Jy. 
A total of 200 galaxies in the Survey satisfy this criteria. 
Most of the selected galaxies are spirals of types Sa-Sc 
while 14 galaxies are mergers, 7 galaxies are close pairs, 
and 3 galaxies have no type determined. In most cases 
we adopt the 60 /xm fluxes from the BGS surveys. For 
galaxies with angular diameters larger than 8', we adopt 
the values from Rice et. al (1988) after multiplying by 
1.18 in order to match the flux scaling (see Devereux & 
Young 1990). This rescaling has little impact on our sam- 
ple completeness limit because most of our sample galaxies 
have flux densities much larger than 5.24 Jy. Some of the 
60 /xm flux density measurements come from Young et al. 
(1989, 2002), and 10 galaxies satisfying our selection cri- 
terion are added to the sample. Five of these galaxies are 
not in the BGS survey area, and two are present in Rice 
et. al (1988) but not in the BGS surveys. 

There are 12 galaxies in in our selected sample that are 
brighter than 5.24 Jy at 60 fim and were not detected in 
CO (see Table 1). The CO (1-0) flux measurements for 
flve of these galaxies are found in the literature. The re- 
maining 7 galaxies are treated as two limiting cases: (a) 
as detections with zero flux; and (b) as detections at the 
upper limit flux value (see Verier 1987). A large number 
of non-detections would severely limit the determination 
of the the true CO LF, but the non-detections account 
for less than 4% in our sample and thus have little overall 
impact. 

Some of the galaxy distances come from direct measure- 
ments such as using Cepheids. We adopt a Hubble con- 
stant of Hq = 75 km s~^ Mpc~^ for the remaining galax- 
ies. Since we are using the CO database obtained from 
largely a single instrument (2.5% of the sample taken from 
other surveys), the internal consistency of the data and the 
analysis should be quite good. Our sample also includes 26 



galaxies in the Virgo cluster. We adopt a uniform distance 
of 16 Mpc for the Virgo galaxies. We have constructed a 
CO LF with and without Virgo galaxies in order to see if 
the presence of the cluster galaxies influenced the shape of 
the CO LF (see § 3). 

To examine the uniformity and completeness of our sam- 
ple, we have analyzed the differential number count statis- 
tics as a function of the 60/um flux density (see Fig. 1). 
If the 60/xm sources are uniformly distributed in a Eu- 
clidean space and arc not evolving, the resulting differen- 
tial number count should be a power-law distribution with 
N oc The number of sample galaxies per bin shown 

in Figure 1 is consistent with such a power-law, but the 
slope is somewhat shallower. When all of the flux density 
bins are included, the best flt power- law index is a = —1.0, 
where log N = No + ax log S. For the 60/Ltm flux density 
between 10 and 100 Jy where most of the sample galax- 
ies fall, the differential number counts arc more consistent 
with the uniform distribution within the statistical uncer- 
tainties. 

There are several factors contributing to the flatter than 
expected power-law slope in the observed number counts 
in Figure 1. First of all, the spatial distribution of galaxies 
is in general not uniform, particularly given the relatively 
high flux density cutoff we adopted. The situation is ex- 
acerbated by the fact that our Galaxy is located within a 
galaxy aggregate called the Local Supercluster. The flat- 
tening of the galaxy counts among the high flux density 
bins as seen in Figure 1 is an immediate outcome of the 
local large scale structures. Soifer et al. (1989) noted a 
similar number count enhancement in their analysis of the 
first Bright Galaxy Sample (BGSl). The mean < V/Vm > 
for our sample is about 0.35, indicating that on average our 
sample galaxies are about 10% closer than the uniform 
case. The < V/Vm > values are particularly low for the 
log LgOfim bins of 9.6-10.4 (see Fig. 2), and similar trend 
is also seen in the < V/Vm > plot by Soifer et al. (1989). 
Comparison to Soifer et al. (1989) value < V/Vm >= 0.47 
is suggesting that major cause of this discrepancy is weak 
bias in our sample towards galaxies closer than average 
expected from the complete sample, while the large scale 
structure plays less important role. Although the FCRAO 
CO Survey selected target galaxies "at random" from the 
parent sample of IRAS and i?-band selected galaxies, a 
slight bias favoring galaxies which were well suited to the 
FCRAO resolution (l'-5' in size) is also present. Such 
a bias manifests non-trivially in these statistics, however. 
Given the limitation of our sample, such as selecting galax- 
ies suitable for observation with FCRAO telescope and the 
local large scale structures, the steep power-law index of 
a ~ —1.2 shown in Figure 1 suggests that our assumption 
of random and uniform sampling from the parent samples 
still seems reasonable. Effects of these potential biases on 
the derivation of the LFs are evaluated by the derivation of 
the 60/im LF using the same sample, as discussed bellow. 



3. LUMINOSITY FUNCTIONS 

A luminosity function represents the space density of 
the sample galaxies per luminosity bin (AL) centered on 
L. It represents the probability density of finding a galaxy 
with a specific luminosity, and it also contains information 
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on the total luminosity in the sampled volume. We use the 
classical 1/Vm method (Schmidt 1968) to derive the LF 
and associated uncertainty directly as: 
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1 



1/2 



(1) 



where p{L) is number of objects per volume per luminos- 
ity bin centered on luminosity L, Vm is the sample volume 
appropriate for each of the galaxies, and N is the number 
of galaxies in the bin. 

The FCRAO Extragalactic CO Survey covers the same 
area as the IRAS BGS (Soifer et al. 1989; Sanders et al. 
1995), and we selected a subset using the same IRAS 60 
/im flux density cutoff. Therefore, the sampling corrccticm 
for our sample can be reduced to a simple scaling relation 
if our sample is a fair subset of the BGS, as argued above 
(§ 2). In this case, the scaling factor is the ratio of the 
number of galaxies contained in BGS surveys and number 
of galaxies in our sample. 

The sample volume obtained using this approach is: 

Vm = 5- = 1-16 X D-" (2) 

where D is maximum distance at which galaxy of given 
luminosity could be detected for a given flux limit, and f2 
is the total solid angle covered by both BGSl (Soifer et 
al. 1989) and BGS2 (Sanders et al. 1995) in sterradian. 
Nfcrao/Nbgs is the ratio of the number of galaxies in 
our sample and the number of galaxies in both BGS sam- 
ples. BGSl and BGS2 together contain 601 galaxies and 
cover 83% of the sky. 

3.1. 60 i^m Luminosity Function 

The 60 iim luminosity represents the contribution from 
the IRAS 60 /xm band to the FIR luminosity (Helou et al. 
1988) and was calculated using the following relation from 
Yun et al. (2001) : 

logLeoi^rniLQ) = 6.014 + 2logD + logSm^Lm (3) 

where D is the distance in Mpc and Seo^im is the IRAS 60 
/Ltm band flux density in Jy. 

As a way to verify our CO sample correction and to eval- 
uate the effects of various sample biases discussed in § 2, 
we first derive the 60 /xm LF for our CO-selected galax- 
ies in order to compare it to the previously published 60 
/xm LF by Yun et al. (2001). As shown on Figure 3, 
the agreement between the two 60 /im luminosity func- 
tions is excellent. Yun et al. sample was selected with 
2.5 times smaller flux density cutoff than the BGS and is 
thus 9 times larger in size. While this is not a complete 
verification of our derivation method, this agreement lends 
support to our assumption that our sample is a fair subset 
of the complete 60/Ltm selected sample. The normalization 
factor could in principle be dependent on the 60;um lumi- 
nosity. The comparison of the 60/Lim LFs suggests that 
this dependence is weak if any. 

About 13% of our sample are Virgo cluster galaxies. 
Again, the 60 /xm luminosity function was constructed in- 
cluding and excluding Virgo galaxies. The resulting lumi- 
nosity functions arc very similar - see Figure 4. The simi- 
larity would arise if the cluster galaxies are distributed in 

^ We adopt the definition Sco = / Icod^, following the same notation as Young et al. (1995). 



luminosity bins in the same manner as the field galaxies 
and if the Virgo galaxies contribute only a small fraction 
to each luminosity bin. Previous studies of Virgo galaxies 
have suggested that their molecular gas contents are af- 
fected little by their cluster environment (e.g. Kcnney & 
Young 1989). Since we are interested in the total cold 
gas density of the local volume regardless of the galaxy 
environment, we include all Virgo galaxies in our analysis 
here on. 



3.2. CO Luminosity Function 

Encouraged by the success of deriving the 60 /im LF 
above for the 200 galaxies in our sample, the CO luminos- 
ity function is derived using the ^^CO(l-O) flux with the 
same sample and volume correction. The CO luminosity is 
calculated as Leo = 4:7rD^Sco, where D is the distance to 
the galaxy in Mpc and Sco^ is the velocity integrated CO 
flux in Jy km s~^. The computed CO luminosity for our 
sample galaxies ranges between 10** and 10* Jy km s~^ 
Mpc^. They are divided into 10 bins of unit magnitude 
(0.4 wide in a log scale). There are larger uncertainties 
at the low luminosity end because only a small number of 
sources occupy these bins. 

The resulting CO luminosity function is shown in Fig- 
ure 4. It is well described by a Schechter function which 
has an exponential cutoff at high luminosity end. The 
low luminosity end power- law slope is cjuitc flat, similar 
to the 60 /im LF (see § 4). To quantify these properties, 
the characteristic parameters are derived for the Schechter 
function (Schechter 1976) of the form 

md{L) = 0*( A)"exp(- A)d(^) (4) 

where (j)* is the normalization factor, a is the low luminos- 
ity end power-law slope, and L* is the characteristic lumi- 
nosity at which the exponential cutoff develops. Since we 
are using logarithmic intervals to construct LF, we need 
to change this equation to an appropriate form: 



p{L) = p*{^r+'expi-^)lnW 



(5) 



in this case p{L) — (j){L)d{L) / d{log{L)) , and it is calcu- 
lated using eq. ( 1). The derived best flt Schechter pa- 
rameters are: p* = (0.00072 ± 0.00035) Mpc'^ mag-\ 
a = (-1.30 ± 0.16), and L* = (1.0 ± 0.2) x 10^ Jy km 
s-i Mpc2 (Af* = 9.4 ± 1.9 X IO^Mq using the standard 
CO-to-H2 conversion factor - see § 4.2). 

The single Schechter function fit using all 10 luminos- 
ity bins produces a rather poor fit with = 50. This 
poor fit is also obvious in Figure 4 as the single Schechter 
function fit shown with a solid line is only marginally con- 
sistent with the data points and their formal uncertain- 
ties. It is possible that a Schechter function is intrinsi- 
cally a poor functional form for the CO LF. On the other 
hand, a Schechter function successfully describes many 
other LFs including the HI LF (e.g. Schneider et. al 1998). 
There is also a theoretical basis in that it is the functional 
form associated with the halo mass function in the Press- 
Schcchtcr formulation of the structure formation (Press & 
Schechter 1974). 
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One possible explanation for such a large value is 
that there are two (or more) distinct populations of galax- 
ies in the sample, similar to the situation with the FIR 
luminosity function as noted by Yun et al. (2001). In 
the FIR, the light-to-mass ratio for the starburst popu- 
lation is systematically enhanced by 1-2 orders of mag- 
nitudes over the field population. The CO luminosity is 
also thought to be elevated significantly among the in- 
tense nuclear starburst systems (Scoville, Yun, & Bryant 
1997; Downcs & Solomon 1998), and this may materi- 
alize as two distinct populations in the CO LF as well. 
There are 14 mergers/starburst systems in our sample, 
and eight of these fall within the last two bins and ac- 
count for about 40% of the galaxies in these bins. When 
a best fit Schechter function is derived using only the 
first 8 bins dominated by normal field galaxies, the 
value drops to 13, which is a significant improvement. 
The Schechter function fit for the bins dominated by nor- 
mal galaxies gives: p* = (0.003 ± 0.001) Mpc'^ mag-\ 
a = (-0.90 ± 0.15), and L* = (3.3 ± 0.6) x 10*^ Jy km s'^ 
Mpc^ (M* = 3.1 ± 0.6 X IQ^Mq). 

The difference between all Schechter parameters for the 
two limiting situations of treating the non-detected galax- 
ies as upper limits or zero detections is smaller than 2%. 
Such a small difference is expected since we have less than 
4% upper limits in our sample. Similarly, including or 
excluding Virgo galaxies has only a minor impact, as dis- 
cussed above (see § 3.1). 

3.3. CO LF from the B-band Selected Galaxies 

As an independent check of our CO LF derived using 

the FIR selection, wc also derive a CO luminosity func- 
tion using an optical -B-band selected sample. The method 
of deriving the CO LF is similar to the method used for 
the infrared selected galaxies: we use the 1/Vm method 
and rescaling of the number of galaxies using a complete 
sample, which is the Revised Shapley-Ames (RSA) catalog 
by Sandage & Tammann (1981) in this case. A declina- 
tion cut off of (5 > —25° and a Galactic longitude limit 
of l&l > 10° are adopted as part of the sample definition. 
This makes the total surveyed area equal to 60% of the 
sky. Our brightness selection criterion is Bt < 12.0™"^. 

Accounting for different Hubble types correctly is a dif- 
ficult issue. Since CO emission is generally very weak or 
undetected among early type galaxies, only the spiral and 
irregular galaxies are included in the analysis. The exam- 
ination of the differential number counts for our B-band 
selected sample yields a power-law index a = —1.35, sug- 
gesting a reasonably good uniformity and completeness. 
The mean < V/Vm > ratio for the sample is 0.33 (see 
Fig. 2), similar to the FIR selected sample. Although we 
adopt the RSA catalog as a complete sample, a few galax- 
ies satisfying our selection criteria are missing in the RSA 
(see Condon 1987). 

The _B-band magnitudes arc taken from the RSA while 
the CO fluxes are taken from the Survey. There are 133 
galaxies in the CO Survey that satisfy the selection cri- 
teria while the parent RSA catalog galaxies satisfying the 
selection criteria are 252 in total. In this case, CO fluxes 
for 13 galaxies (10%) are upper limits. The CO luminosity 
function derived for this i?-band selected sample is shown 
on the Figure 5. This LF is very similar to CO LF from 



the FIR selected galaxies but with a lower characteristic 
luminosity L* and a larger scatter on the low luminosity 
end. The lower characteristic luminosity L* is a direct 
consequence of the absence of mergers in this sample. The 
best fit Schechter parameters for the i3-band derived CO 
luminosity function are: p* = (0.0021 ± 0.0009) Mpc'^ 
mag-i, a = (-1.0 ± 0.2), and L* = (4.8 ± 1.1) x 10*^ Jy 
km s-i Mpc2 {M* = 4.5 ± 1.0 x IO^Mq). 

4. DISCUSSION 
4.1. Uncertainties and biases in the CO LF 

One major source of uncertainty in the derived CO lumi- 
nosity function is that it is derived indirectly using the FIR 
and optical B-band selection functions. A direct deriva- 
tion from the observed CO properties is possible in princi- 
ple, but the complication associated with the a priori un- 
known line widths adds a significant uncertainty. Instead 
wc take advantage of the known tight correlation between 
FIR and CO luminosity in deriving the CO LF for our 
sample. The 1.4 GHz radio luminosity function derived 
by Yun et al. (2001) using the radio-FIR correlation and 
IRAS 60 yum flux density agrees very well with the radio 
LFs derived directly (Condon, Anderson, & Helou 1991; 
Condon et al. 2002), giving some assurance to this tech- 
nique. We have examined whether any systematic trends 
(thus a bias) exist in the Seoi^m/Sco ratio as a function of 
the 60 pm flux density and luminosity. As shown in Fig- 
ure 6, little trend is seen in this ratio, and the known linear 
correlation seems to hold well within the uncertainties for 
the entire range of flux density and luminosity. 

Deriving a luminosity function for one wavelength using 
a selection function at another wavelength may be robust 
enough to work even if the two quantities arc correlated 
in a non-linear way. There is a known correlation be- 
tween FIR, CO, and optical B-band luminosity for late 
type galaxies, but the correlations involving the B-band 
luminosity are shown to be non-linear (sec Young et al. 
1989; Perea et al. 1997). Yet, the CO luminosity func- 
tion derived using the FIR selection function (Figure 4) 
is in excellent agreement with the CO LF derived using 
the B-band selection function (Figure 5). The selection 
bias is present as some of the most CO luminous galax- 
ies (merger starbursts) are missed in the B-band selected 
sample. Nevertheless the agreement is striking given the 
non-linear dependence between CO and B-band luminos- 
ity. 

Another source of a significant uncertainty is the unifor- 
mity of sampling and correctly accounting for the selection 
function for the sample of galaxies which is not complete 
by the adopted selection criteria. Both the differential 
source count and the < V/Vm > analysis suggest a non- 
uniform distribution of the sample galaxies (see § 2), and 
a slight bias towards galaxies brighter at 60pm may be re- 
sponsible for this effect. The large scale structure also have 
influence in lowering < V/Vm > value for our sample. The 
non-trivial nature of these effects are demonstrated by the 
fact that including and excluding Virgo cluster galaxies 
from the sample make little difference to the differential 
source count and the < V/Vm > values. 

4.2. Local Cold Gas Mass Density 
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CO emission is a commonly used tracer of molecular hy- 
drogen because CO is one of the most abundant molecules 
in cold ISM and because its excitation in astrophysi- 
cal conditions is determined by collisions with hydrogen 
molecules. Using the CO LF derived above, we can for the 
first time estimate the total mass density of molecular gas 
in the local volume. Adopting N{H2)/I{C0) = 3 x 10^° 
cm~^ [K km s""'^]""'^ (see review by Young & Scoville 1991), 

M{H2) = 1.18 X IQ'^ScoD^Mq (6) 

where Sco is total CO flux of a galaxy in Jy km s~^ and 
D is luminosity distance in Mpc. 

The molecular mass density in the local volume con- 
tributed by each luminosity bin is shown in Figure 7. The 
dominant contribution to the mass density comes from 
galaxies around L* as expected. The summation over the 
10 bins gives Ph^ = Y^MHjVm = (2.4 ± 0.7) x IO^Mq 
Mpc~^. Integration of the LF using the Schechter param- 
eters obtained in previous section gives a little bit smaller 
values: pn^ = (2.2±1.1) x lO'^Mg Mpc^^ ^j. ^^i^ g^- trough 
ah 10 bins, and pn^ = (2.2 ± 0.9) x IO^Mq Mpc'^ for fit 
trough first 8 bins. The second fit is more realistic, since 
it fits much better bins that dominate contribution to the 
total mass, i.e. bins around L* . The uncertainty stated 
is ±1(7. The systematic uncertainty, which includes un- 
certainties in the flux measurements, distance determina- 
tions, and CO-to-H2 conversion, is probably larger. Unless 
otherwise is stated, we adopt pH^ = (2.3 ± 0.9) X lO^M© 
Mpc^'^ as an average value between values obtained from 
the fit and the direct summation. Since the dependence of 
the gas mass density on the Hubble constant {h = Hq/IOO 
[km s~^ Mpc~^]~^) is linear, this result can be written as 
Ph^ = (3.1 ± 1.2) X IO^/iMq Mpc--\ 

The molecular hydrogen mass density in the local vol- 
ume obtained from the i?-band selected sample is ph2 = 
(3.1 ±0.9) X 10^/iMq Mpc-3 using a direct summation of 
the contribution of each galaxy and (3.1 ± 1.5) x IO^/zMq 
Mpc~^ using the best fit Schechter parameters. Values for 
the local molecular mass density obtained from the FIR 
selected sample and the B-band selected sample are in 
good agreement. Since we have better statistics for the 
larger FIR selected sample, we adopt the local molecular 
gas mass density derived from the FIR selected sample 
here on. Using observed II2/HI mass ratio for different 
morphological types of galaxies and fraction of each mor- 
phological type Fukugita et al. (1998) estimated a similar 
value for pH2- 

In comparison, Zwaan et al. (1997) estimated the local 
atomic gas mass density of pHi = (5.8 ±1.2) x 10^ HMq 
Mpc~^ from their Arecibo HI strip survey while Rao & 
Briggs (1993) derived pHi = (4.8±1.1) x IO'^HMq Mpc'^ 
from a sample of optically selected galaxies. Therefore, 
the molecular gas mass density in the local volume is about 
50-65% of the atomic mass density, and molecular gas rep- 
resents a significant component of the total mass density 
of the neutral gas. The HI masses for 176 galaxies in our 
sample are known (Young et al. 1995), and we have com- 
puted the HI mass density in the local volume using the 
same procedure as for the H2 mass density. A summation 
over the HI mass weighted by Vm gives a value consistent 
with the value obtain by Zwaan et al. (1997). 

We examined our sample for any trends in Mh^/Mhi 
ratio vs. Mhi and Mh^- For the range of HI masses, 



8.8 < log Mhi(Mq) < 11, the ratio of molecular to neu- 
tral atomic hydrogen masses is on average around unity 
with considerable scatter (see Fig. 8a). Masked in the large 
scatter is a possible trend in M^^/Afe/ as a function of 
Hubble type (see Yoimg et al. 1989). A clearer trend of 
increased Mh^/Mhi ratio with increasing H2 mass is seen 
in Figure 8b. For log MH^i^o) < 8.5, the average of this 
ratio is below 0.5 while the ratio jumps to around 2 near 
Mh.2 ~ IO^Mq. Above H2 mass of lO^M©, on average 
galaxies have more moleciilar than atomic gas. 

Using the molecular gas mass density derived here and 
the value of HI mass density from Zwaan et al. (1997) we 
derive a total cold gas mass density at present epoch as 
a fraction of the critical density VIH1+H2 = (3-2 ± 0.8) x 
lQ~'^h~^. Supposing that He contributes 25% of the total 
gas mass density, we add 33% to the derived value of ilu, 
which gives 9.HI+H2+He = ^gas = (4.3 ± 1.1) X 10~^/i~^. 
These values for total neutral gas mass density are 15% 
smaller if we use Rao & Briggs (1993) value for the HI 
mass density. 

The present baryon density estimated from the D/H ra- 
tio and Big Bang nucleosynthesis is fif, = 0.02 ± 0.002/i~^ 
(Buries et. al 2001) while the estimate from the cos- 
mic microwave background anisotropy is also around fib ~ 
0.02/i"^ (de Bernardis et. al 2002). The stellar mass con- 
tents account for fi* ~ 0.0025ft.-^ (Cole et al. 2001; but 
also see Benson, Frenk, & Sharpies 2002). We conclude 
that the cold gas content of late type galaxies is around 
20% of the stellar mass content and about 2% of the total 
baryonic content in the current epoch. 

4.3. CO-to-H2 conversion 

The lack of electric dipole moment makes direct observa- 
tions of molecular hydrogen difficult in general, and study- 
ing the spatial extent and molecular gas mass requires an- 
other tracer. Highly abundant, chemically robust, and 
easily excited by collision with H2 molecules, CO is the 
most commonly used tracer of molecular gas. The CO (1- 
0) transition is optically thick under most astrophysically 
interesting conditions, making it relatively insensitive to 
metallicity and abundance effects. The two key excitation 
parameters of density and temperature have a nearly can- 
celing effect, making CO a fairly reliable tracer of H2 in a 
broad range of physical conditions (see reviews by Maloney 
& Black 1988; Young & Scoville 1991). 

The derived ratios of N{H2)/I{CO) range between (1 — 
5) X 10^0 cm-2 [K km s''^]''^ (Bloemen et al. 1986; Dick- 
man et al. 1986; Scoville & Sanders 1987). We adopt a 
constant CO-to-H2 conversion factor of N{H2)/ I{CO) = 
3 X 10^° cm~^ [K km s~^]~^ (sec discussions by Young 
& Scoville 1991). Devereux & Young (1990) show that 
the H2 mass estimates of galaxies from their CO luminos- 
ity are accurate to ±30%. Young & Scoville (1991) show 
that the CO-to-H2 conversion for galaxies of diverse mor- 
phology and metallicity are similar in absolute value to 
the conversion in the Milky Way. While the H2 mass esti- 
mate for an individual galaxy may be uncertain to about 
30%, the H2 mass estimate for an ensemble of galaxies 
should be more reliable. In most galaxies, giant molecular 
clouds (GMCs) and cloud complexes dominate the total 
molecular gas mass, and adopting a conversion factor con- 
sistent with the values derived from the Galactic GMCs 
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should yield the most robust, mass- weighted estimates of 
gas masses. 

Among the low metallicity galaxies such as Small Mag- 
ellanic Cloud, CO abundance may become low enough to 
affect the standard assumption of self-shielding and ther- 
malization. The derived conversion factors are generally 
larger (e.g. Wilson 1995), and low metallicity dwarf galax- 
ies are often undetected entirely in CO. Therefore the H2 
mass and mass density derived from the CO luminosity are 
strictly lower limits since molecular gas from these galax- 
ies is missing from our analysis. On the other hand, the 
contribution by the low luminosity bins to our derived H2 
mass density is small (sec Fig. 7) as the majority of the 
contribution to the total mass density comes from bins 
around L* . Therefore, we can largely neglect the effect of 
low metallicity, low luminosity galaxies in the derivation of 
the total molecular gas mass density for the local volume. 

As discussed already in § 3.2, the CO emission may be 
elevated for the FIR luminous galaxies with nuclear gas 
concentrations, and the CO-to-H2 conversion factor may 
be smaller than the canonical value. In the extreme en- 
vironment of the nuclear starburst regions, CO emission 
arises from a multi-phase medium with sub-thermal exci- 
tation in the diffuse phase, and the standard conversion 
factor may over-estimate the molecular gas mass by a fac- 
tor as large as 3 to 5 (Scoville, Yun, & Bryant 1997; 
Downes & Solomon 1998). Such merger/starburst galax- 
ies contribute significantly only for the top two luminosity 
bins in our CO LF, and these two bins contribute less than 
5% to the the total molecular gas mass density. Therefore, 
the high CO luminosity conversion factor docs not have 
a significant impact on the derived total molecular mass 
density. 

Similarly, increased CO emission in the central 1 kpc of 
our Galaxy and other galaxies has been suggested by re- 
cent observations (Paglione et. al 2001, S. Hiittcmeister, 
private communication) . High angular resolution observa- 
tions of gas-rich spiral galaxies frequently reveal a distinct 
component in the central 1 kpc (see Sakamoto et al. 1999), 
and some fraction of the total CO luminosity may arise 
from such a component. On the other hand, among the 
Virgo spirals studied by Kenney & Young (1988), a simple 
exponential distribution without a significant central com- 
ponent offers a good fit for the 12 out of 14 galaxies whose 
CO emission is spatially well resolved (i.e., CO detected at 
> 4 positions along the major axis). We made no attempt 
to account for enhanced CO emission in the central kpc as 
the required information is generally not available. If the 
molecular gas properties of Virgo spirals are typical of the 
late type field galaxies, then the possible contribution by 
the enhanced nuclear CO emission in some galax;ies may 
not be substantial. 

In summary, using the canonical CO-to-H2 conversion 
relation is problematic in some cases, such as low metal- 
licity systems or luminous nuclear starburst systems, and 
CO is a poor tracer of molecular gas among low metallic- 
ity dwarfs. For these reasons the local molecular gas mass 
density we infer is really a lower limit. However, since the 
majority of the mass contribution comes from L* galax- 
ies whose molecular gas mass is dominated by GMCs like 
our Galaxy, the use of the standard conversion factor still 
offers a fairly reliable estimate of the total molecular gas 



density for the local volume. 

5. SUMMARY 

Utilizing the largest available CO survey, the FCRAO 
Extragalactic CO Survey (Young et al. 1995), the CO 
luminosity function for the local volume is derived us- 
ing (1) a FIR selected sample of 200 galaxies that satisfy 
S&Ofim > 5.24J?/ and (2) optical B-band selected sample 
of 133 galaxies. Although neither of the samples is com- 
plete in terms of the sample selection, a sampling function 
is constructed using a well defined parent sample, and a 
non-parametric CO luminosity function is derived from 
each sample. By examining the properties of the CO lu- 
minosity functions, we conclude: 

1. The CO luminosity functions derived from the FIR 
and B-band selected samples are reasonably well 
described by a Schechter function. The character- 
istic luminosity L* is around Leo ^ 10^ Jy km s^-'^ 
Mpc^. The low luminosity end of the CO luminos- 
ity function is (a = —1.3 to — 0.9). Similar values 
are obtained for the two CO LFs derived using the 
two differently selected samples. 

2. The molecular gas mass density of the local volume 
is = (3.1±1.1) X 10^/iMq Mpc"-'' which is about 
50-65% of the HI gas mass density. This value is 
not strongly affected by variations in the CO-to-H2 
conversion factor for low and high luminosity galax- 
ies since it is dominated by the L* galaxies with 
Leo ^ 10' Jy km s'^ Mpc^ {M^^ - 5 x IO^Mq). 
Its dependence on the global CO-to-H2 conversion 
factor is linear and should be secure to within a 
factor of 30% or better. 

3. Combined with the HI gas mass density, we esti- 
mate the total cold gas mass density at the present 
epoch as a fraction of the critical density Q,hi+H2 = 
(3.2 ± 0.8) X lO-^/j-^ When the He contribution 
is included, the cold gas mass density increase to 
^Hi+H.+He - ^gas = (4.3± 1 . 1) X IQ-^/i" 1 . There- 
fore, the cold gas content of late type galaxies cor- 
responds to about 20% of the stellar mass content 
and about 2% of the total baryonic content in the 
universe. 
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Fig. 1. — Number of galaxies per bin (0.2 in log scale) of 60 /im flux density and integrated number of galaxies. Although the number of 
galaxies per bin grows with decreasing flux density, the power-law slope of a ~ —1 suggests a non-uniform distribution of sources. The dotted 
line represents the expected behavior of uniformly distributed non-evolving sources in Euclidean geometry, appropriate for the 601 galaxies in 
the two BGS samples. The fact that high flux bins axe higher than the complete sample line is noticed also in BGS (Soifer et al. 1989) and 
is caused by large scale structures. The lowest luminosity bins are lower than expected probably because sources near the flux cutoff that are 
spatially resolved by IRAS have systematically lower flux entries in the Point Source Catalog and are thus missed. 
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Fig. 2.— V/Vm analysis for (a) the 60 /im LF and (b) the CO LF. The mean value of V/Vr 
on average about 10% closer than the uniform case. 
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Fig. 3. — IRAS 60 /xm LF obtained using our CO sample of 200 galaxies corrected for the sample selection. The compajrison with the 60 
/im LF derived from a much larger sample obtained by Yun et al. (2001) shows a good agreement. 
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Fig. 4. — The CO luminosity function, derived using the same sample and volume correction as the IRAS 60 fim LF shown in Figure 3. The 
solid line is the best fit Schechter function derived using all 10 bins. The CO LF follows a Schechter function shape closely with a faint end 
power-law index a = —1.30 and the characteristic luminosity L* = 9.8 X 10^ Jy km s~^ Mpc^. The dotted line shows the best fit Schechter 
function derived excluding the two highest luminosity bins (see § 3.2). The power-law index for this fit is a = —0.90 and L* = 3.3 x 10® Jy 
km s~i Mpc^. 
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Fig. 5. — The CO luminosity function obtained for the optical _B-band selected sample derived from the RSA catalog as the parent sample. 
The solid line shows the best fit Schechter function (see text). 
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Fig. 6. — (a) Average S^o^rn/Sco vs. LeOMm- High luminosity bins have higher values of SeOMm/ScOi but they also suffer from the small 
number statistics, (b) Average Seo/xm/Sco vs. Seo/xm- The absence of any clear trend in either plots suggests that our IRAS selection does 
not introduce any obvious bias to the CO measurements. 
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Fig. 7. — Molecular mass density contribution per magnitude of H2 mass. The H2 mass bins near the M* value, log ~ 9-5, dominate 
the contribution to the total mass density. The solid and dotted lines correspond to the same lines shown in Figure 4. 
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Fig. 8. — The molecular to atomic hydrogen mass ratio are shown as a function of (a) Mhi and (b) Mh2- Individual galaxies are shown as 
triangles while the bin averages are shown with filled circles and connected in solid lines to accentuate the trend of growing M{H2)/M(HI) 
ratio for larger H2 masses. 
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Table 1 

Galaxies with Sso > 5.24 Jy not detected in the Survey 



Name 


R.A. 

[h m s] 


dec 
[d m s] 


D 

[Mpc] 


[Jy] 


Sco 
[Jy km/s] 


Type 


Ref. 


NGC 337 


00 59 50.3 


-07 34 44 


23.7 


9.33 


140 


Sd 


u 


NGC 470 


01 19 44.8 


+03 24 35 


35.6 


7.09 


180 


Sb 


u 


NGC 693 


01 50 31.0 


+06 08 42 


21.9 


6.86 


70 


so/a 


u 


NGC 925 


02 27 17.0 


+33 34 43 


9.5 


9.03 


566 


Sd 


1 


NGC 3227 


10 23 30.6 


+19 51 54 


20.6 


8.32 


756 


SAB (pec) 


* 


NGC 3486 


11 00 23.9 


+28 58 30 


7.4 


6.24 


263 


Sc 


1 


NGC 4027 


11 59 28.4 


-19 19 55 


19.5 


11.89 


60 


Sdm 


u 


NGC 4526 


12 34 03.1 


+07 41 59 


16.0 


5.63 


87 


so 


2 


NGC 4532 


12 34 19.3 


+06 28 07 


16.0 


8.93 


60 


Im 


u 


NGC 4808 


12 55 48.9 


+04 18 13 


16.0 


6.92 


100 


Scd 


u 


NGC 6824 


19 43 40.9 


+56 06 33 


49.0 


5.94 


30 


Sb 


u 


Mrk 273 


13 44 47.5 


+55 54 11 


153.8 


22.09 


62 


Merg. 


3 



Coordinates are for J2000 and are taken from NASA Extragalactic Database, ^eo is IRAS 60 /um flux density taken from 
BGSl and BGS2 except for NGC 925 where we used value from Rice et. al (1988) increased by 18%. References for Sco 
(column 8) are: l.-Sage (1993); 2.-Sage & Wrobel (1989); 3.-Sanders et. al (1991); u - upper limit from Young et al. 
(1995). *-a new measurement of central 45" obtained using the FCRAO 14-m telescope in March 2002. NGC 3227 was 
not included in the Survey originally because this galaxy was confused with NGC 3226. 



